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Abstract. We apply our quantum macrostatistical treatment of irreversible processes 
[18] to prove that, in nonequilibrium steady states, (a) the hydrodynamical observables 
execute a generalised Onsager-Machlup [12] process and (b) the spatial correlations of 
these observables are generically of long range. The key assumptions behind these results 
are a nonequilibrium version of Onsager regression hypothesis [11], together with certain 
hypotheses of chaoticity and local equilibrium for hydrodynamical fluctuations. 
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1. Introduction 

It is now well appreciated that a key problem in the statistical mechanics of irreversible 
processes is the characterisation of nonequilibrium steady states [2,6,13], and a number of 
different approaches have been made to this problem. Some have employed rather general 
strategies, based, for example, on a hypothesis of Anosov dynamics [6,13]; while others 
comprise treatments of concrete microscopic stochastic dynamical models [20,4,1]. 

A different approach to nonequilibrium statistical mechanics has been made by the 
present author in a number of works [16-18] that are centred on the hydrodynamical ob- 
servables of quantum systems. This approach, like Onsager's [11] treatment of the subject, 
is designed to form a bridge beween the microscopic and macroscopic pictures of mat- 
ter, rather than a deduction of the latter from the former. Its basic assumptions concern 
only very general, model-independent properties of many-particle systems, and its scope 
is thus intended to be complementary to that of works based on microscopic treatments 
of many-body problems. The results to which it has led [18] include a mathematical 
characterisation of local ther mo dynamical equilibrium and a generalisation of Onsager's 
reciprocity relations to a regime where the macroscopic dynamics is nonlinear. 

In this note, we extend the macrostatistical scheme by the introducation of a chaoticity 
hypothesis for the fluctuations of the non-conserved currents associated with the locally 
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conserved hydrodynamical observables in nonequilibrium steady states. On this basis we 
obtain the results that, in these states, 

(a) the fluctuations of the hydrodynamical observables execute a generalised Onsager- 
Machlup (OM) process [12]; and 

(b) the spatial correlations of these observables are generically of long range. 

The latter result constitutes a mathematical generalisation of results previously proved 
for certain special classical stochastic models [20,4,1]. At a heuristic level, similar results 
concerning long range correlations have also been obtained from Landau's fluctuating hy- 
drodynamics [8,5]. 

We remark here that the result (b) marks a qualitative difference between equilibrium 
and nonequilibrium steady states, since the hydrodynamical correlations of the former 
states are generically of short range*, except at critical points. 

2. The Model 

The Quantum Picture. We take the model to be an A^-particle quantum system, 
S, that occupies an open, bounded, connected region, Q^, of a d-dimensional Euclidean 
space X and is coupled at its surface to an array, 7?., of reservoirs. We assume that the 
particle number density, v, of S is AT- independent and that fijv is the dilation by a factor 
Ljy of a fixed, A^-independent region ft, of unit volume. Thus 0,^ = L]yQ={LNx\xEfl} 
and Lj\r = (N/v)^/'^. In a standard way, we represent the observables and states of E by 
the self-adjoint operators and density matrices, respectively, in a separable Hilbert space 
Hn- We assume that, as has been established under rather general conditions [14, 22], the 
composite system (E -|- TZ) evolves to a unique steady state, a;jv, as t— >oo. We denote the 
expectation value, for this state, of an observable A by {ujn'iA). We shall assume that all 
interactions are invariant under space translations and rotations. 

We assume that E has a finite set of linearly independent, extensive, conserved observ- 
ables Q = (Qi, ■ -jQn)? which intercommute up to surface effects and are thermodynami- 
cally complete [18] in the sense that the states corresponding to pure equilibrium phases 
are labelled by the expectation values of their global densities in the limit N^oo. We 
denote by s{q) the equilibrium entropy density corresponding to the value q of the global 
density of Q in this limit. 

We assume the observables Q have locally conserved, position dependent densities 
(^qi{x), . ., qnix)) := q{x), and we denote their cvolutes at time t, in the Heisenberg picture 

of the dynamics of the composite (E -|- TZ), by q^^\x,t)=q^^\x) = [q[^^\x), . .,q^^{x)). 
In accordance with the standard requirements of quantum field theories [21], we assume 

that the fields qj^^ are operator- valued distributions. Thus, denoting by V{Qn) the space 
of continuous linear functionals on the subset of the L. Schwartz space 'D{X) [15] with 

* Our distinction between 'long' and 'short' range will be expressed in a sharp mathe- 
matical form in Section 4 
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support in Qn, we assume that the fields q^^^ are operator-valued elements of 'D'{Qn)^- 

The Hydrodynamical Picture of E. We assume that this is given by a continuum 
mechanical law governing the evolution of a set of locally conserved classical fields qt{x) = 
.,qn,tix)), on a macroscopic space-time scale that we shall presently specify. 
Here the fields qj,t{x) represent the densities at position x and time t of the extensive 
thermodynamical variables of S. In general, their dynamics is given by a classical equation 
of motion of the form 

the right hand side representing the value of a functional JF of qt at the position x. For 
simplicity we base our treatment here on the case where this equation corresponds to a 
nonlinear diffusion and thus takes the form 

^ = V.{K{q,)Vqt), (1) 

subject to certain stationary, spatially varying boundary conditions determined by the 
reservoirs TZ. Here and elsewhere qt is a function of position, and K{qt) is an n-by-n 
matrix [Kjk{qt)], whose action on Vqt is given by standard matrix multiplication. We 
shall presently specify the relationship between qt and the quantum field qt^\ 

We assume that the unit of the macroscopic length scale is the distance Ln, introduced 
above; and correspondingly, since Eq. (1) is invariant under the scale transformations 
we take the macroscopic time scale to be L^. Thus, Eq. (1) represents a 
phenomenolgical dynamics in the bounded spatial region O on this time scale. We assume 
that that equation admits a unique stationary solution, qt — q, subject to the imposed 
boundary conditions. We assume throughout that the system is in a single phase region 
wherein both K{q) and the entropy density s{q) are smooth functions of q as this variable 
runs over the range of q{x) for x&Q. 

A simple important consequence of Eq. (1) is that, although the fields qt are locally 
conserved, the associated currents are not. This is crucial for our key assumption, formu- 
lated in Section 3, to the effect that the irreversibility of the flow of qt stems from chaoticity 
properties of the fluctuations of these currents*. 

Relationship between the Quantum and Hydrodynamical Pictures. We assume that the 
classical field q is the expectation value of q for the steady state ujn-, as represented on the 
macroscopic length scale of unit Ljv. Thus, 

q{x) = limjv^oo('^Ar; qiL^x)) V xeQ. (2) 

Correspondingly, since the unit of the macroscopic time scale is L%, we represent the time- 
dependent fiuctuations of qt{x) (=q{x,t)) about its steady state value, on the macroscopic 



* In the case of Navier-Stokes hydrodynamics, it would be the stress tensor, not the 
mass current, that would be the locally non-conserved field, and our chaoticity assumption 
would pertain to the fiuctuations in this tensor. 
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scale, by the field 

= N'/^[q(''\L,,x,L%t) - {cON;q{LNx))], (3) 
the normalisation parameter N^/'^ being the canonical one for fluctuations [7]. 

3. The Fluctuations: a Generalised Onsager-Machlup Process. 

We now assume, on grounds [16-18] representing the intercommutativity of the quan- 
tum fields qj,t{x) and qk,u{y) at macroscopic space-time separation and the fact that the 
ratio of macroscopic to microscopic relaxation times becomes infinite in the limit N^oo, 
that the quantum process ^^^^ reduces to a classical stationary Markov process in this 
limit, i.e. that 

Vxi,. .,a:^eO, ti,. .,tm^Ii+, (4) 

where ® denotes the R"^ tensor product. Since V spaces are complete, it follows from this 
formula that is a distribution- valued random field. The forward time derivative of ^t, 
as defined by Nelson [10], is 

D^t = limh^+oh-^E{^t+h - 616), (5) 

where E{.\^t) denotes the conditional expectation functional, given the field 6- 

Regression Hypothesis. We now invoke a version of Onsager's regression hypothesis 
[11], to the effect that the regressions of the fiuctuations 6 ^ire governed by the same 
dynamical law as the 'small' perturbations 5qt of the macroscopic field q. To this end, we 
start by inferring from Eq. (1) that the linearized equation of motion for Sqt is 

g-/qtix) = C5qt{x) := W .{K{q{x))W{5qt{x)) + [K'{q{x))Sqt{x)]Vq{x)), (6) 

where K'{q) is the derivative of K{q), i.e. its gradient w.r.t. the variable q: thus 
[K'{q)6qt]jk — Yl7=i[^^jk{q)/dqi]Sqi^f We assume that the perturbation of q does not 
change its boundary conditions, i.e. that Sqt vanishes on the boundary of Q. 

We assume that Sqt, like the quantum field q, is a D'(fi) "-class distribution and that 
the operator C, defined in Eq. (6), is the generator of a one-parameter semigroup of 
transformations {Tt|teR+} of as determined by the equation 

-^=m;T„=I; (7) 

and thus that the solution of Eq. (6) is 

5qt = Tt-sSqs V t>s>0. (8) 
4 



We assume that the system satisfies the dissipativity condition that the perturbation 5qt 
vanishes in the hmit t— >oo, i.e. that 

hm^^ooTt^ = 0, (9) 

for all elements of that vanish on the boundary of Q,. 

We now assume, as a generalisation of Onsager's regression hypothesis for equilibrium 
fluctuations that, for t>s>0, the conditional expectation of ^t, given ^g, is 

E{it\is)=Tt-sL\ft>s>Q. (10) 

Hence, by Eqs. (5), (7) and (10), 

D^t = C^f (11) 

Further, by Eq. (10) and the stationarity of the process, the two-point function 
E{^t{x)®^t'{x')) takes the value {Tt-t'®I)E{^{x)®^{x')) or {1®%' -t)E{i{x)®^{x')) ac- 
cording to whether or not t>t' . Thus, 

E{it{x)®^t'{x')) = 

{Tt-t'®I)E{ax)®ax'))0{t - t') + iim'-t)E{CixMix')) (1 - eit - t')), (12) 

where 9 is the Heaviside function that takes the value 1 or according to whether or not 
its argument is non- negative. 

Extended Stochastic Process: the Currents. Recalling that the fields qt are locally 
conserved observables, we assume that the above formulation of the stochastic process 
has a canonical extension to a larger process (Ct^Vt), where % represents the fluctuations 
of the currents* associated with qt. Thus, and % conform to the local conservation law 

^ + ^-Vt = 0, (13) 
and consequently, by Eqs. (5) and (11), 

^t = -E{V.rit+o\^t). (14) 

We now deflne 

fit = r]t- E{r]t+o\^t) (15) 

* Note that it follows from the local conservation laws that the divergences of these 
currents are also classical. Hence, it is irrelevant for our purposes that the currents them- 
selves might not be classical, since they enter into our calculations only through V.ry. In 
fact, we have obtained the same results by a fully quantum treatment [19] of the currents, 
which does not contain formal conditional expectations such as that in Eq. (15). 
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and 

d 

h = -^■fjt= - (16) 

where and r/^^t are the ^u'th Cartesian components of x and rjt, respectively. We note 
that Eq. (15) may naturaUy be interpreted as signifying that E(r]t+o\^t) and f]t are the 
secular and residual stochastic parts, respectively, of and thus Eq. (16) signifies that 
bt is minus the divergence of the stochastic part of the current. 



It follows now from Eqs. (13)- (16) that 

= ^t + bt, (17) 



dt 

which has the form of a Langevin equation*. Our interpretation of it, however, will depend 
on the chaoticity and local equilibrium hypotheses that we shall presently introduce. At 
all events, bt is statistically independent of $, (:= ^o) for t > 0, since, by Eqs. (7), (10) and 
(17), 

mio = ^t^m) - mm = (| - = o. 

Further, by Eq. (17), 

E{bt{x)®bt,{x')) = C®I){— ~ I®C')E{^t{x)®^t,{x')), 

where C is the image of £, as defined by Eq. (6), under the transformation x^x' . Hence, 
by Eqs. (7) and (12), together with the identity de{t - t')/dt=S{t - t'), 

E{bt{x)®bt,{x')) = -[E{Ci{x)®i{x'))+E[^{x)®C'i{x'))]5{t-t'). (18) 

The higher order correlation functions of bt will be governed by the following chaoticity 
hypothesis. 

Chaoticity Hypothesis. This is the hypothesis that the space-time correlations of the 
stochastic part of the non-conserved current fiuctuations associated with as viewed on 
the microscopic scale, are of short range. Since the ratios of the macroscopic to microscopic 
scales of both length and time are infinite, this signifies that the space-time correlations of 
the currents rit{x) have zero range. Further, since the fiuctuations of quantum fields with 
short range correlations are generally Gaussian in the large scale limit [7] , we assume that 
fjt is a Gaussian process. Thus, our chaoticity assumption is that fjt is a Gaussian field 



* Strictly speaking, fit and bt should be treated as distributions w.r.t. t, since it will 
presently emerge that each of them corresponds to a white noise (cf. Eqs. (18) and (29)). 
Mathematical propriety can easily be achieved, however, without changing the structure 
of the argument, by working with the integral wt = Jq dubu, which would correspond to a 
Wiener process, instead of bt. 
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with zero range space-time correlations. Hence, by a standard theorem on distributions 
[15; Theorem 35], we can state the chaoticity hypothesis as follows. 

( C) fit is a stationary Gaussian process whose two- point function E {fji^^t{x)<^fj^,t' {x')) 
is a finite linear combination of S{x — x')S{t — t') and its derivatives, with coefficients given 
by generalised functions of x. 

We note here that the Gaussian property of bt implies that of ^t, for the following 
reasons. Since C is the generator of the semigroup T, it follows from Eq.(17) that 

pt+to 

Ct+to = Tt+toC + / dsTt+to-sbs V t, to>0. 
Jo 

By (C), the integral on the r.h.s. of this equation is Gaussian, while, by Eq. (9), the first 
term on its r.h.s. vanishes in the limit t^oo. Hence, as the stationarity of the process 
implies that it is isomorphic with Ct+to^ it must be Gaussian. 

Equilibrium Conditions. As a preliminary to the formulation of local equilibrium 
conditions for the nonequilibrium steady state, we first formulate the true equilibrium two- 
point functions of the fields and For this purpose we assume that the equilibrium 
state is achieved by arranging the reservoirs TZ so that q{x) is uniform on the boundary, 
dQ, of Q. With this boundary condition, Eq. (1) has a stationary solution in which q{x) 
is uniform throughout O. Consequently, by our uniqueness assumption of Section 2, q{x) 
reduces to a constant, q in the equilibrium state. Further, recalling our assumption of 
Section 2 that the interactions of the system are translationally and rotationally invariant, 
we assume that the corresponding symmetries are not broken in the pure equilibrium phase 
and thus that the process (^t, fjt) is invariant under the space translations and rotations 
that are implementable within the confines of O. We remark here that the limitation in 
the Euclidean symmetry imposed by the boundedness of O is not serious from the physical 
standpoint, since any point of this open region corresponds, in the microscopic picture, to 
one that is infinitely far from the boundary of the quantum system S. 

Under the above assumptions, a quantum statistical treatment of the field ^ yields 
the following result for the static two- point function for ^ (cf. [18, Ch. 7, Appendix 
C]): in fact, it corresponds to a thermodynamic limiting version of the Einstein formula 
P = const. exp(5'). 

EequiliCixMix')) = J{<1)5{X - X'), (19) 

where E^quii is the equilibrium expectation functional and J{q) is minus the inverse matrix 
of the Hessian of the entropy density function i.e. 

J{q) = -[d^s{q)/dqjdqk]-\ (20) 

On the other hand, it follows from the assumptions of translational and rotational invari- 
ance that the two-point function of 'f]t{x) takes the form 

Eequil{fl^i,t{x)®flT,,t'{x')) = G{X -x',t- t')6,^^, (21) 
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where G is a scalar w.r.t. the vector space X. Hence, the chaoticity condition (C) reduces 
here to the following form. 

{C)equil- G{x — x',t — t') IS a finite linear combination, with constant coefficients, of 
the distribution 5{x — x')S{t — t') and its derivatives. 

In order to pin down the explicit form of G we now note that, by Eq. (6), the constancy 
of q implies that C is simply K{q)A, and therefore, by Eqs. (16), (18) and (19), that 

Eequil{bt{x)<^bt'{x')) = Q^f Eequil{fji^,t{x)<^fi^,t'{x')) 

= 2[Kiq)Jiq)lyrr,ASix - x')5{t - t'), (22) 

where, for any matrix A, [Ajs^^ denotes the arithmetic mean of A and its transpose*. 
Hence, by Eqs. (21) and (22), 

AG{x -x',t- t') = 2[K{q)J{q)]sym^S{x - x')5{t - t'). 

This equation, together with condition {C)equih fixes the form of G according to the 
formula 

G{x -x',t- t') = 2[K{q)J{q)]symS{x - x')5{t - t'), 
and consequently, by Eq. (21), 

Eequil {rinA^)®riv,t' {x')) = 2[K{q)J{q)]symS{x - x')5{t - t')5^^. (23) 

We now note that, for any positive e and any xqGO and to^R'+j Eqs. (19) and (23) 
are invariant under the transformations x-^xq + ex, x'^xq + ex', ^— ^^e*^/^^ and x^xq + 
ex, x'^xq + ex' , t— >to + e^^, t'^to + e^t' , fjt^e^'^^/'^f]t, respectively. Thus they are 
equivalent to the following formulae. 

e'^Eeguiii^ixo + ex)®C{xo + ex')) = J{q)S{x - x'), (24) 

and 

e'^'^^Eequii{fjfj.,to+eH{xo + ex)<^fii,^to+eH'{xo + ex')) = 2[K (q) J{q)]sym5{x - x')5{t-t') . (25) 

Evidently, if e is chosen to be 'small', then these equations correspond to local conditions 
concentrated at the space-time point {xo,to). 

* In fact, it has been proved in [18] that K{q)J{q) is a symmetric matrix, i.e. that the 
Onsager reciprocity relations prevail, subject to the assumption of microscopic reversibil- 
ity. Moreover, this result was extended there to the nonequilibrium situation under the 
assumption of a certain local equilibrium hypothesis not employed in this note. However, 
we shall not assume this symmetry here, since it is not needed for our present purposes. 
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Local Equilibrium Conditions. This last observation leads us to propose the following 
local equilibrium conditions for the general stationary nonequililibrium situation. 

lim,^oe'^E{^{xo + ex)®C{xo + ex')) = j{q{xo))5{x - x') (26) 

and 

lime^oe'^^'^ E{fi^^to+eH{xo + ex)^f]^,^to+eH'{xo + ex')) = 

2[K{q{xo))j{q{xo))]^^J{x - x')5it - t')6^,. (27) 

The Generalised Onsager-Machlup Process. In view of the chaoticity condition (C), 
it follows immediately from this last equation that the two-point function for -qt must take 
the following form, since the presence of derivatives of S{x — x')S{t — t') would render the 
l.h.s. of Eq. (27) divergent. 

E{v^4x)(S)V.,Ax'))^2[K{q{x))j{q{x))]^^J{x-x')S{t-t')5^^ (28) 

It follows from this formula and Eq. (16), the two-point function for bt is 

E{bt{x)^btix')) = -2V.{[K{q{x))j{q{x))]^^^V5{x-x'))5{t-t'). (29) 

Hence, by the Gaussian assumption in (C), bt is a white noise and consequently, by Eq. 
(17), executes a generalised Onsager-Machlup process [12]. We remark here that the 
two-point function for bt given by Eq. (29) is of the same form as that of the noise term 
in Landau's fluctuating hydrodynamics [9]. 

4. Long Range Spatial Correlations. 

The static two-point function for the fluctuation field ^ is 

W{x,x') = E{^{x)^^{x'))=E{^t{x)®^t{x')) by stationarity. (30) 

Since the ratio of the macroscopic to microscopic length scale is infinite, short range cor- 
relations on the latter scale reduce to zero on the former one. Accordingly, we term the 
range of the correlations 'short' or 'long' according to whether or not it reduces to zero 
in the macroscopic picture. Hence our condition for long range correlations is simply that 
the support of the distribution W does not lie in the domain {{x,x')&fl'^\x = x'}. The 
following Proposition establishes that the spatial correlations of ^ are generically of long 
range. 

Proposition. Let $ be the n-by-n matrix- valued functional on the classical field q 
defined by the formula 

^q; x) = A [K{q{x))j{q{x))] + V.^{q; x)sym (31), 
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where 



[Jmk{q{x))^qi{x) - Jik{q{x))Wqm{x)]. (32) 

Then under the above assumptions, a sufficient condition for the spatial correlations of ^ 
to be of long range is that $(g) does not vanish. 

Comments. (1) The Proposition establishes that the correlations are generically 
of long range, since the condition that ^{q) vanishes can be satified only for special re- 
lationships between the functions K(^q{x)) and By contrast, the corresponding 
correlations for equilibrium states are generically of short range, except at critical points. 
A treatment of critical equilibrium correlations of fluctuation observables is provided in 
Ref. [3]. 

(2) In the particular case of the symmetric exclusion process [20,4], n = 1, d = 
1, K{q) = 1, s{q) = —qlnq — (1 — g)ln(l — q) and q{x) = a + h.x, where a and h (t^O) 
are constants. Thus, in this case, it follows from Eqs. (20), (31) and (32) that $(q';x) = 
— 26^^0. Hence, as proved by other methods in Refs. [20, 4], long range correlations 
prevail in this model. 

Proof of Proposition. Suppose that the spatial correlations of ^ are not of long 
range, i.e. that the support of W does lie in the domain {{x.,x')&Q?'\x — x'}. Then it 
follows from a classic theorem on distributions [15; Theorem 35] that W{x, x') is a finite 
linear combination of 5{x — x') and its derivatives, with coefficients given by generalised 
functions of x. Under this assumption, it follows from Eqs. (26) and (30), together with 
the continuity of J(g(a;)), that this combination must reduce to the form 

W{x,x') = J{q{x))6{x-x'), (33) 

since the presence of derivatives of 5{x — x') would cause the l.h.s. of Eq. (26) to diverge. 
On comparing Eqs. (18) and (29), we see that 

{C®I + I®C')E{^{x)mx')) .{[K{q{x))j{q{x))]^^^V5{x - x')). (34) 

Hence, as a®h is the transpose of 6(8)a, 



*jfc(g;a;) = 



o 



Kji{q{x)) 



[{mi)E{C{xM{x'))],^^ + [{C'^l)E{C{x'M{x))]^^^ = 

2V.{[K{q{x))j{q{x))]^^^V5ix - x')), (35) 
i.e., by Eqs. (30) and (33), 

[C{j{q{x))Six - x'))]^^^ + [C'{j{q{x'))Six' - x))] = 

2V.{[K{q{x))J{q{x))]VSix - x')). (36) 
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We now infer from the formula for C, given by Eq. (6), together with Eq. (32), that 

A[K{q{x))j{q{x))S{x-x')] +V.[^{q;x)5{x-x')], (37) 
while the interchange of x and x' in this formula yields the equation 

C'{j{q{x'))5{x' - x)) = A'[K{q{x'))j{q{x'))5{x' - x)] + V'. [^(q; x')S{x' - x)] = 

A' [K{q{x))j{q{x))5{x - x')] + V. [*(g; x)5{x - x')] = 

K{q{x))j{q{x))A5{x' -x) - ^{q; x).V5{x - x'). (38) 

It now follows easily from the last two equations, together with the definition (31), that 
the difference between the left and right hand sides of Eq. (36) is simply $(g; x)S{x — x'). 
This signifies that the assumption of short range correlations implies that ^{q) = 0; or 
equivalently that a sufficient conditions for the static ^-correlations to be of long range is 
that $(g) does not vanish. 

5. Concluding Remarks 

The quantum macrostatistical theory presented here is based on the regression hypoth- 
esis together with the assumptions of local equilibrium and chaotic current fluctuations. 
On this physical basis, we have obtained both a generalised Onsager-Machlup process and 
a generic picture of long range correlations in nonequilibrium steady states for systems 
whose phenomenological dynamics corresponds to a multi-component nonlinear diffusion. 
Our derivation of these results depended on a coordination of the spatial and temporal 
macroscopic scalings, which was facilitated by the invariance of the assumed phenomeno- 
logical law, given by Eq. (1), under the scale transformation x^Xx, t^XH. We remark 
that, for systems whose phenomenological laws do not have any simple scale invariance, 
e.g. for Navier-Stokes hydrodynamics, the situation is clearly much more complex; and a 
generalisation of our results to such cases would presumably require an intricate multi-scale 
analysis. 
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